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INTRODUCTION 

In  this  paper  I  summarize  some  results  obtained  in  the  period 
June  198S  -  July  1986.  It  is  my  pleasure  to  acknowledge  collabora¬ 
tion  with  the  following  colleagues:  M.J.  Ablowitz,  P.  Clarkson,  M. 
Kruskal,  R. A.  Leo,  L.  Martina,  U.  Mugan,  V.  Papageorgiou,  P.M.  Santini, 
and  G.  Sol  iani . 

“tU  The  results  on  Inverse  Scattering  in  multidimensions  and  on  the 
algebraic  properties  of  equations  in  2+1  (i.e.  two  spatial  and  one 
temporal)  dimensions  should  be  of  particular  interest:  With  respect  to 
algebraic  properties  of  equations  in  2+1  wo  note  that,  the  question  of 
finding  the  recursion  operator  and  the  ^-Hamiltonian  formulation  of 
these  equations  has  remained  open  for  a  rather  long  time.  It  was  even 
doubted  in  the  literature  if  the  relevant  results  in  1+1  could  be  ex- 
tended  to  2+1.  Santini  and  the  author  have  recently  shown  that 

equations  in  2+1  solvable  via  the  Inverse  Scattering  Transform  are  bi- 
Hami 1  ton lan  systems.  T bey— also  given  m p 1 i c i t the  recursion 
and  bi -Hami 1 tonian  operators  for  large  classes  of  equations  in  2+1, 
including  the  Kadomtsev-Petviashvi 1 i  (a  two  dimensional  analogue  of 
the  Korteweg-deVries)  and  the  Oavey-Stewartson  (a  two  dimensional 
analogue  of  the  nonlinear  Schrodinger)  equations.  _ 

In  this  paper  I  emphasize  the  basic  ideas  and  results.  Further¬ 
more  an  attempt  is  made  to  put  these  results  into  perspective.  De¬ 
tails  can  be  found  in  the  sited  papers. 
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1.  AN  INVERSE  PROBLEM  FOR  N  x  N  AKNS  IN  MULTIDIMENSIONS 

This  problem  has  been  studied  in^,  where  I  considered  the  in¬ 
verse  problem  associated  with  the  following  system  of  N  first-order 

equations  in  n+1  dimensions: 
n 


f  +o  l 
x0  1=1 


qf. 


a  =  oR  +  ioj,  Oj  t  0,  n  >  1,  (1.1) 

where  q(xQ,x)  is  an  N  x  N  matrix-valued  off-diagonal  function  in  Rn+*, 
decaying  suitably  fast  for  large  x^.x,  and  the  are  constant  real 

diagonal  N  x  N  matrices  (we  denote  the  diagonal  entries  of  J.  by 

IN  * 

Alternatively,  using  the  transformation 

n 

T  (x0,x,k)  =  u  (x0,x,k)exp[i  £  k|l(x)l-  oxQJ^  )],  k  c  I  ,  (1.2) 

I  considered  n  1 

u w  »0£  (J,u  +  ikJj-.u])  «  qu.  (1.3) 

x0  1*1  *■  xi  *■  1  1 

I  assume  that  n  <  N,  otherwise  the  entries  of  the  matrices  will  be 

linear’y  related  and  one  can  always  reduce  n  by  a  change  of  coordin¬ 
ates.  An  inverse  problem  in  this  case  is  defined  as  follows:  Given 
appropriate  inverse  data  T,  where  T  is  an  N  x  N  matrix-valued  off- 
diagonal  function  of  suitable  inverse  parameters,  reconstruct  the  po¬ 
tential  q. 

There  is  a  twofold  motivation  for  considering  such  an  inverse 
problem. 

(a)  If  o  *  -1  then  the  above  reduces  to  the  formulation  of  a 
physically  important  inverse  scattering  problem:  Given  the  scattering 
amplitude  S ( A , k ) ,  A,  k  c  Rn,  which  is  a  function  of  the  scattering 
parameters  A,  k,  reconstruct  q. 

(b)  In  recent  years  a  deep  connection  has  been  discovered  bet¬ 
ween  inverse  scattering  of  linear  eigenvalue  problems  in  one  spatial 
dimension  and  the  initial  value  problem  of  certain  nonlinear  evolution 
equations  in  1+1  (i.e.,  one  spatial  and  one  temporal  dimension).  Re¬ 
cently  a  similar  connection  has  been  used  to  extend  the  above  results 

to  nonlinear  evolution  equations  in  2+1  (i.e.,  two  spatial  and  one 
2-6l 

temporal  dimension)  .  In  particular  the  inverse  scattering  of  (1.3) 
with  o  *  -1  and  n=l  has  been  used  to  linearize  the  N-wave  interaction 


equations  in  2+1  (see  Ref.),  the  Davey-Stewartson  (DS)  I  (see  Ref.  J) 
(a  2+1  analog  of  the  nonlinear  Schrodinger) ,  and  the  modified  Kadomt- 
sev-Petviashvi 1 i  (MKP)  I  (see  Ref.4^)  (a  2*1  analog  of  the  modified 
KdV)  equations.  Furthermore  the  inverse  problem  of  (1.3)  with  o  =  i  and 
n*l  has  been  used  to  linearize^  DSII  and  MKPII.  However,  in  spite 
of  the  above  success  in  2+1,  no  physically  interesting  equation  is 


known  to  be  related  to  (1.3)  for  n  >  1  and  a,  t  0  [the  N-wave  mterac- 

71 

tlon  equations  in  n+1  spatial  and  one  temporal  dimension  J  are  related 


to  (1.3)  but  with  Oj  =0], 

The  novelty  associated  with  inverse  problems  in  greater  than  two 
spatial  dimensions  (n  >  1)  stems  from  the  fact  that  while  the  potent¬ 
ial  q(xg,x)  depends  on  n  +  1  variables,  the  inverse  data  Tfkp.kj.m^, 


.  ,m  ),  k„  e  Rn,  k,  c  Rn,  m  c  R,  depends  on  3n  -  1  variables.  This 
n  k  l  c 


has  important  implications: 

(a)  The  inverse  data  must  be  appropriately  constrained.  This 
“characterization"  of  the  inverse  data  is  conceptually  analogous  to 


the  characteri zation  of  the  inverse  scattering  data  in  the  multidim- 

8- 1 1 1 

ensional  Schrodinger  equation 


(b)  The  existence  of  “redundant"  scattering  parameters  in  the 

inverse  scattering  of  the  Schrodinger  equation  is  used  to  reconstruct 

the  potential  in  closed  form  in  terms  of  the  scattering  amplitude 

12  1 

function.  This  is  the  well-known  Born  approximation  J.  Can  one  use 
the  redundancy  of  the  inverse  parameters  here  to  also  reconstruct  q 
in  closed  form? 

In  the  above  paper,  I  do  the  following. 

(a)  Following  A.  Nachman  and  M.J.  Ablowitz,  I  derive  an  equation 
that  characterizes  inverse  data:  d  , .  d  , .  n!jp[T](w  w..P') 

-1  *  T'J/-  w  ,1  -if  _ 5 _ ! - - - - - .  M  41 


COPY 

k INSPECTED 

X  i 


d'i  n!^[t](w  w. 


(wQ,w)  %  T1J(w0,w,<) 


where  wQ  cR.  w  c  Rn,  <  <Ln'1  are  related  to  k.  m  and  N  is  a  quadrat- 


fUAstff 


ic  function  of  T.  That  is,  T T  -*  {  k ,  m )  is  appropriate  inverse  data  iff 
the  right-hand  side  of  (1.4)  is  independent  of  *.  Hence,  Eq.  (1.4)  serves 


as  both  characterizing  T1^  and  defining  t'j. 


(b)  1  reduce  the  general  problem  of  recons  true 1 1 ng  an  N  x  N  po- 


^n/ _ 

tty  Codes 
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tential  q  in  n  ♦  1  dimensions  to  one  of  reconstructing  a  2x2  poten¬ 
tial  with  entries  q1J,  qJ1  in  two  dimensions.  The  inverse  data  needed 
for  this  reconstruction  is  precisely  T1^,  T31.  This  reduction  makes 
crucial  use  of  the  existence  of  redundant  scattering  parameters.  In 
this  sense  it  is  the  analog  of  the  Born  approximation.  However,  the 
crucial  difference  is  that  while  in  the  inverse  scattering  of  the 
multidimensional  Schrodinger  equation  one  can  reconstruct  the  potential 
in  closed  form,  here  one  can  only  reduce  the  general  problem  to  one 
for  2x2  matrices  in  two  dimensions.  This  reduced  problem  was  solved 

in4l 

Partial  results  about  the  case  a  *  i  were  given  in13^.  Equation 
(1.3)  was  also  considered  in14^  where,  although  the  characterization 
problem  was  solved  (an  equation  very  similar  to  (1.4),  the  so-called  “T 
equation",  was  first  obtained  in34^)  the  problem  of  an  effective  re- 
constructon  of  q  was  left  open. 

The  basic  steps  are  as  follows: 

1.1  Bounded  Eigenfunctions. 

The  function  u(xg,x,k)  defined  below,  solves  equation  (1.3),  is 
bounded  for  all  complex  values  of  k  and  tends  to  I  for  large  k: 

(Ojj|)r  exp[iB1J(x0-^0,x1-E1,k)] 

‘  „2 


I 


u  J(xg,x,k)  =  6 


’J 


♦sgn- 


2*i  j  • 


(xl^l) 


°Ji(V5o) 


*  (qu)ij(E0.S1.x2-(x1-E1)j’/ji . xn_(xl'f;i)Jn/Jl  *k>  ’  k  e  (1’5) 

where  81J  is  defined  by 


lj'*  ,Xj ,k)  *  2  [r  '-,2l 

1  1*1  °I  l 

es 

sgn(ojjj) 


B  J(x 


Equivalently  u i j  satisfies 


X1  (° k  ^)  I 


k  *k  Hk  .(1.6) 
R  I 


u1J(x0,x,k)  *  JRn+l  d^od^cn-l]Rn-l 

exp[i8i3(x0-S0,x1-£1,k)](qu)13(E0,  E.k) 


dm1 


Zgia^x-E.m)^ 


KreroJ!(Veo) 


(1.7) 


where 


dm  *  dm,. ..dm.  a  (x.m)  i  Z  m  (x.  -  x.  ~r) ,  c  i  - (1.8) 

2  n  i«2  1  1  1  jJ  n  (?Tr)n 

1.2  Departure  from  Hoi omorphi city. 

Let  u,J  be  defined  by  eq.  (1.5).  Then 

— (xn.x,k)  *  Z  y’u!  -  ji)exp[1B,j(xn.x1.k)] 


P  P 


*  cn-lj  n-1  ^exp[ia1(x,m)]TlJ(k.m)u(x0.x.A1J(k,m))Eij.  (1.9) 

R 

where  81 J(xn,x. ,k) ,  a’fx.m)  are  defined  by  (1.6),  (1.8)  respectively; 

E.  .  is  an  N  x  N  matrix  with  zeros  in  all  its  entries  except  the  ij  , 

J  i  i  i  i 

which  equals  1;  and  A  J  and  T  J  are  given  by 

,i 

l]j(k,m)  i  (kjJ  -  Z  m  4,  k.  ),  ^j(k,m)  =  (k  +m  ,k  );  r=2,..n. 

1  1r  4=2  Ji  1  R  1 

y1  ♦o/4'ri!jJo1|. 

Tlj(k,m)  *  [  n+1  dC0d^expC-i  S1  J  (C  0 .  C  j ,  k )  -  i  ct1  ( C  ,m )  ]  ( q  u) 1 J  ( CQ .  fi.  k ) .  ( 1 . 10) 


1.3  Characterization  of  T. 

(a)  Assume  that  3 u/3 Sp  is  given  by  Eq.  (1.9)  and  the  T1J(k,m) 
is  given  by  (1.10).  Then 

l’yj(k,m)  ■  •  Icn  „  .dM2TU(/j(k,M),m-K)T*j(k,m) 

*  -  (jJ-jj)(Jp-Jp)]  *  N’j[T](k.m).  (1-11) 


.U  *  (J‘  .  jj)  _L  -  (j1  .  jJ)-L  . 

rp  P  P  r  r  3E 


(1.12) 


(b)  Assume  that  is  given  by  Eq.  (1.9)  and  that  a  u/dRrdli 

r  j  j  r 

is  symmetric  with  respect  to  r,p.  Then  T  J(k,m)  solves  (1.11). 

Following  A.  Nachman  and  M.J.  Ablowitz  1  introduce  appropriate 
Born  variables.  Then  equation  (1.11)  can  be  integrated.  Furthermore, 


we  can  compute  the  limit  of  T  J  in  the  new  coordinates  as  | Xp |  -  ® 
( see  below) : 

Let  wjj.wjj.wr  £  »  2 . n  e  R1  and  c  (E1,  £  =  2,...,n,  be 

defined  by 


n  J 1  -J J  .  n  J1 -JJ  n  J1 

"0J  *  Z  ~Ll^\2k-r  .  wjJ  *  -  E  — p(ok  )j  -  £n  -j  , 
U  r- 1  J  \  rI  1  r=l  Ojjj  r  1  r=2  r  jJ 


w£  •  m £’  *V  •  Tj-TT*  4  =  2 . n- 

JrJi 


(1.13) 


Assume  that 


(Jj-J^)(Jp-Jjj)  +  ( Ji“Ji ) ( Jp“Jp) *  for  a11  distinct  i.j.r  Ahd  p (*1. (1.14) 


for  convenience  of  writing  we  usually  suppress  the  superscripts,  i , j 
in  Wg.Wj,  x.  Let  k  denote  kj .... ,kn,  m  denote  m^, ... ,mn,  x  denote 
*2.....Xn,  w  denote  w^,..,w  .  Then  we  have  the  following. 

(a)  The  inverse  of  the  transformation  k,m  -  wQ,w,x  is  given  by 

k^X^i-oj), 1*2, kj  =  -  MJ^-jj.)xr  ♦ 


(a / 1  a i  2,-  *  c"  w  j’ 

+ -  o  -r-i  u:  .  a. 15) 

Ji  -  Ji 

(b)  In  the  new  coordinates,  Eq.  (1.11)  with  r*l  becomes 
jt' 3  i  i 

— — (w0,w,x)  *  Njj[T](w0,w,x),  p  *  2 . n.  (1.16) 


(c)  In  the  new  coordinates 

Tlj(wn,w,x)  *  (  p+1  dE0dEexp[-i(wnEn+wE)](qu)lj(En,  E,wn,w,x), 


0’  '""O' 


where  <  !  I  vf  . 

r=l  r  r 


(1.17) 


(d)  Let 


“  H«  u*j.  (1.18 

1  U  U  I  |  ,  i  1 


Then  the  uf3  satisfy 


,2n 


ui  (*o,x,wO’w'  =  sgn  2tTT  cn-l 

dx^dx 'dm2exp[i {(x0-Xg)w0+w0+(x-x' )w}]  . ii#  ,  . 

- -■ -  q1J(xn,x  )^J'|J(xn,x,  ,wn,w), 

Vi  -  «J  (v*o>  0,00 


.ii  sgn(o.J') )  , 

ui  (x0*x*w0*w)  *  1  *  STi  cn-l  J 


,2n 


(1. 


dx’dx’dm  q-' 1  (xi,x‘  )jj^(xi,x’  ,wn,w)  . 

— - —< - 3 - - - - - ,  G*J=0,  for  all  i,  Iti,  If  j. 

xrxi  ‘  aJl<x0"x0) 


(e)  lim  Tlj(w0,wlX)  =>  j  n+1  d£0d£exp[-i  (w0£0+we;)  ] 


*  q1J(t0,OC^(£0.S.w0,w)  *  Tlj(w0,w). 


(1. 


(f)  The  basic  characterization  equation  is  given  by 

lf  dvvllJp[T1(v",xP  ) 


Tl;,(w0,w)*T1J(w0,wtx)  -  ,jrj  2 


XP  ’  XP 


(1.21 


where  xP  denotes  x2 . Xp-i ,xp,xp+l ' 


1.4  Reconstruction  of  q. 


It  follows  from  the  above  that  as  1  X p  I  ■*'00 »  the  u1J,s  decouple. 
Furthermore,  the  u]J,  u'j'1  satisfy  a  system  of  two  equations  depend¬ 
ing  on  q1J,  qJ1.  It  turns  out  that:  a)  By  introducing  appropriate 


spatial  variables  £,  the  u3J,  u'j'3  satisfy  equations  in  two  spatial 


equations,  b)  The  invers  data  needed  to  reconstruct  u!J,  u^J  (and 
hence  q1J,  qJ1)  can  be  obtained  from  Tig  : 


Let 


.J,1  ,J,i 

a  ♦  SJLlL-1 


jIjJ-jJj1 


1  r  ulur 


J1J2*J1J2 


r  i1  iJ  iJi1 
J1J2'J1 J2 


r  *  1 , . . . ,n. 


(1.22 


where  for  convenience  of  writing  we  have  suppressed  the  dependence  of 
ar.  8  on  i , j .  let  tn  £  R,  L  c  Rn, 


Then  we  have  the  following. 

(a)  The  system  (1.19)  becomes 

*  s9nT7T  f 


,2  d«odcif*i-*i  ‘  oJ!(co  -  s^r1 


exp[iBlj(C0-^),';1  -  q.kllq^a^UQ.q,^ 


-  . q.k), 


■  1+sq"  °-fy  L  d«odqCcrq  -  oJ^Q-^)]*1 


qJ1hJ(C0lCi>C2  •  (€r€i)J2/Jl*53 . Vk,» 


(1.24) 


~  n  ^2*^2  i 

k  t  t  (k  a  ♦  4-4  k  fl  ),  81J(x0,x.,k) 

r=  1  r  r  jJ-jj  r  r  U  1 

-ij/  -  2  (ok). 

8  (x0.*j.k)  »  4_J  [x0|o|2kj  -  xj  -j4]. 

1 

(b)  T,J  in  the  new  coordinates  becomes 


(1.25) 


T1J(k.m)  =  n+1d^dCexp[-iB1J(^.q,k)  ♦ 


+  "»2(52  '  Si  "71  +  E 

q  r=3 


(1-26) 


m2  *  m2  +  z  mrsr,  =  ny,  l  =  3,...,n. 


(1.27) 


(c)  The  inverse  data  associated  with  (1.24)  and  the  analogous  problem 
for  ui’.ui’  are  given  by  f’^t  TJ ' .  Let 


T,J(k.52-qjJ/j;.C3 . C„)  *  cn-l)Rn-l  dm  mpC'W4!  JT} 


+  i  I  m  £  ]T1J(i,fn). 
r=  3 


(1.28) 


f’Ju.q-qJ^/jj.q . £.n)  *  2  dC'dqexp[-i61J(q,q.k)] 


(q^K^.q.q  -  (q-q)jJ/jJ,q,...,cn,k). 


(1.29) 


2.  INVERSE  SCATTERING  FOR  THE  HYPERBOLIC  N  x  N  AKNS  IN  MULTIDIMEN- 


This  problem  has  been  studied  in15J,  where  I  considered  equations 
(1.1)  and  (1.3)  with  o  *  -1.  This  system  appears  to  be  physically 
more  interesting  than  (1.1)-(1.3):  (a)  Since  the  system  is  hyperbolic 
one  may  consider  the  physically  important  question  of  inverse  scat¬ 
tering  (IS);  i.e.,  given  a  scattering  amplitude  function  S (A , k )  find 
the  potential  q(xQ,x).  (b)  A  special  case  of  the  above  system,  namely 
if  the  Jjs  are  constrained  by 
j  ^ 

— ?,  p,r*l,...,n,  i . J =  l . N,  (2.1) 

-r  jJ-jJ 

i  is  associated  with  the  N-wave  interaction  in  n*l  spatial  and  one  tern- 

*•*'  y  igl 

poral  dimensions  ’  . 

With  respect  to  (a),  (b)  above  the  following  results  are  obtain- 

ed: 

(a)  I  first  define  S(A,k),  »,  k  e  Rn,  in  terms  of  an  eigenfunc¬ 
tion  uL(xQ,x,k) .  S( A ,k)  motivates  the  introduction  of  the  Born 
variables  wQ  c  R1,  w  e  Rn,  x  e  F  .  I  then  define  T(A,k)  in  terms 
of  an  eigenfunction  u(xg,x,k).  The  crucial  new  step  is  that  u  is  re¬ 
quired  to  have  analyticity  in  one  of  the  x’s.  5ay  x^*  as  opposed  to 
one  of  the  k's.  Let 


54 


then  t'j  satisfies 


j  f’j(Vw)  =  Tij(wQ,w,y)-(P[  T1  J  )  (wg.w., )  , 

where  P  denotes  a  (+)  or  (-)  projection  in  the  variable 
'l 

(p;  f)(x2)  :=  !  {  d>-f(x^)/2ri(x-  -  x2  ;  iO). 

2  ^  JR 

The  sign  of  certain  parameters  1 ,  where 


(2.2) 
i .  e . , 


e£rj1  :=  (J1  -  ‘  Jr)(Jl  *  J1J*  (2-3) 
determines  whether  the  (+)  or  the  (-)  projection  is  needed.  Equation 
(2.2)  defines  T1J,  which  actually  depends  on  q1J,  and 


:iJ  :=,  u* 

I X2 I -« 

The  question  whether  (2.2)  is  sufficient  for  the  character! zation  of 
T1J  remains  open.  With  a  proper  coordinate  transformation,  u33, 
-j1,  define  an  IS  problem  for  a  2  x  2  potential  with  entries  q1J, 

qJ1,  in  two  spatial  dimensions.  The  inverse  data  needed  for  this 

problem  are  simply  related  to  T1J,  TJ1.  The  solution  of  such  an  IS 
41 

problem  was  given  in  J.  It  is  interesting  that  there  exists  a  simple 
relationship  between  T  and  S.  Actually  if  N  *  3  then  T  =  S.  This  is 
quite  remarkable  since  for  the  first  time  a  closed-form  expression 
can  be  obtained  between  S  and  T.  Furthermore,  u  L  satisfies  a  Volterra 
as  opposed  to  a  Fredholm  equation;  hence  one  ironediately  excludes  the 
possibility  of  bound  states. 

With  a  proper  coordinate  transformation,  the  N-wave  interaction 
equations  in  n  ♦  1  spatial  dimensions  can  always  be  reduced  to  two 
spatial  dimensions.  Thus  a  genuine  three-spatial -dimensional  nonlin¬ 
ear  evolution  equation,  related  to  an  IS  problem,  remains  to  be  found. 

Let  uL(xg,x,k),  Xg  e  R^,  x,k  e  Fn  be  the  solution  of  (1.3)  which 
also  solves 

-L^(x0,x,k)  =  *  f  0  dx^exp[ i k ( J  -J3  ) ( xQ-Xg) ] 

;  -oo 


(quL)lJ(xQ,x  +  J1(xq"Xq) ,k) , 


where  — 


t  k  J‘,  mx  =  Z  m  x  c  :«  (2-t )  “  , 
rs  1  r=l 


i  ?  3 

and  x  +  J  xQ  denotes  Xj  +  JjXg,...,xn  +  Jnxo'  Let  *(xq.x)  be  the  gen¬ 
eral  solution  of  (2.1)  with  o  =  -1  tending  to  f(x  +  XqJ)  or  g(x+XQj) 
as  Xq  -  -<»  or  Xq  -  respectively.  The  scatteirng  operator  is  de¬ 
fined  by  g  =  Sf,  and  S  is  uniquely  defined  in  terms  of  S(X,k),  where 

s'^(X.k)  :=  cnj  n+j  dXgdx  exp[i(k-x)x  +  i (kJ^-XJ ’ )xQ] 

»  (q  ) 1  j(x0,x,k),  i  t  j,  (2.5) 


kx  :  =  Z  k  x.  . 
i«l  *'  '• 

The  Born  variables  are  defined  by 

j  i  kz 

wQ  kJJ  -  XJ  ,  wt  :«  kt  -  *  4.  i  »  2 . n,  (2.6) 

Jf  Jl 

where,  for  convenience  of  writing,  we  suppress  the  dependence  of  wQ, 
v  on  i ,  j . 

T 1 J  ( * ,  k )  satisfies  an  equation  similar  to  S1-1  where  is  re¬ 
placed  by  u|  .fXThe  eigenfunction  satisfies  an  equation  similar^to 
(2.4)  where  |  ®  is  replaced  by  !  .  This  integral  is  either  j  ^ 

"*  oil 

or  -  j  according  to  the  following  requirements :  (1).  If  -  0 

‘O  'Xq 

then  choose  j  .  (2)  If  f  0  then  (xq-Xq)  must  l1dve  the 

same  sign  for  all  l  (i,j  are  fixed).  (3)  If  (2)  can  be  satisfied 

'*0  fxo 

with  either  ;  or  -i  ,  choose  the  first  integral.  To  illustrate 


the  above,  consider  N  *  4.  Since 


,  -  .  fji  „ 


s  0,  there  exist  only 


two  nonzero  ^  ,  if  they  are  of  the  same  sign  choose  both  integrals 

*  *  0  '  *  o 

to  be  ,  if  they  are  of  opposite  signs  then  one  integral  is  j 


and  the  other  is  - 1 


In  the  Born  variables,  T.^(wq,w,<)  depends  on 


u-r 
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^J{xQ,x,w0.w,x  ) ,  which  satisfy 


-Ji 


X  Q  ,  X  ,  Wq  ,  W  ,  X ) 


iji 


dx^expt 1 [w0-wJ 


r-2 


c;J\r] 


(xn-x')}(qu.)<-J(x'  x  +  Jl(xn-x' ),wn,w,<). 


0  0 


0  O'’  0* 


(2.7) 


Equation  (2.7)  implies  the  following:  (1) 


ij 


and  hence  T .  .  have 


l’  1  0i 

analyticity  properties  with  respect  to  X2:  Y  J(o^J)  is  a  (  +  )  or  (-) 
function  with  respect  to  X?  (i.e.,  it  is  holomorphic  in  the  upper  or 
lower  x2  half-planes,  respectively)  according  to  whether  (xQ-Xg)  e  £ 
is  >0,  or  <0.  (2)  As  i  x2  |  c  -  u- J  (w0,x,w0,w) ,  where 

uS>  satisfy  a  reduced  system:  -  0  for  all  i  t  l,j,  and  u  ^ 

satisfy  a  system  of  two  integral  equations  of  the  Volterra  type  (see 
below) .  Hence  as  i  v  2 1  -  ”, 

T,J-  f,J(w0,w)  :«  cn[  j  dxQdx  exp[i'w0x0  ♦  wx)] 

4  K 


'j 


xn,x)u'JJ(xn,x,wn,w) . 


Since  T1^  is  a  (  +  )  or  (-)  function  of  x2  tending  to  T1J  as  lx2l-  ®. 
its  (-)  or  (  +  )  projection  must  satisfy  Eq.  (2.2).  [We  define  p  (l)*j, 
p'(l)  -  - 

Given  T1^,  T^1  one  can  compute  T1^,  T^  1  ,  which  actually  can  be 
used  to  reconstruct  the  2x2  matrix  potential  with  entries  q''1,  q^  : 
We  consider  the  reduced  system.  The  crucial  fact  is  that  it  corre¬ 
sponds  to  N  =  2,  and  hence 

J2Jr~JrJ2  *  k(JJ-J1), 

jJ  -  J 


JJr=  or  J’  ar 


k  :*  ^  -V-;  (2.8) 

,1  ,J  ,J  ,1  ,J  ,1 

J1J2'J1J2  J1  ‘  J1 

1 ,  j  *—  i,  and  again  for  convenience  of 
writing  we  suppress  the  ij  dependence  of  ar>  i  r-  Since  one  can  intro¬ 
duce  a  single  k,  it  follows  that  the  reduced  system  must  be  transform¬ 
able  to  a  system  in  two  dimensions.  This  is  indeed  the  case.  Let 


6  is  defined  as  a  with  2 
r  r 


v  s r  x=f  <  :f  -v*l;  x  = 

0  ’0*  1  V  2  ’2  Y  t>  “  n*  a 
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'  (?2->1jC0  ♦  51J5i).  i  -  3 . n.  | 

where  ylj  :■  (j}jd  -  jjjJ)/(j{  -  j{).  i’j  :»  (J2  -  jj)/(j[  -  j{). 


Then  the  reduced  system  satisfies 

,  x 

*  6^  +  cnJ  0  d''  i  .  d r’ dm  exp:  i [ k ( J f-J1? )  +  m.jM  « 

,  n  ...  <2-10> 

«  Uo-^H^jjUrq)  +  ir:2n,rur-5r,i(q-i),tJ(56*5,*k)' i=i>j> 

where  (qui ) 1  ^  =  q1^^,  (q£.)3j  =  qJ1"’3.  Hence  2^  satisfy  two  inte¬ 
gral  equations  in  the  variables  Lg,  \  ^ : 


n  ^3  ( r  r  N)  s  i'-J+c  '  ^  j.  i  | 

“ i  U0”l’  '  -  clj  °^OjR 


f  Aq  ^ 

*clj  U  d't)i  2  d^idmiexP'  1  [k{ J )  +  nijjj 


[[(S0-^)  +  im1Jj-(t,1-q)  •  ( q  u  i )  J(^g,q,i,2 . ^n,k). 


(2.11) 


The  inverse  data  associated  with  (2.11)  can  be  obtained  from  T1J:  Let 


w  j=w  .,  i  =  3 , . . .  ,n ,  w2=  w 


Then,  since 


WQXg  *  WX 


2"  w2  *  :,Vr-  '=  Z,l'r*  2r}V  (2'I2) 

r=j  r=l 


(kjd  -  ■  Jj ) --0  +  (k-  -Kj  *  2  w  •;  , 

r  =  2 


T  in  the  new  coordinates  becomes 

Tld(k,X,w)  =  cn  i  n+1  dCgdfexpf  i  (  kjd  -  AjJjfg  +  i  ( Ic- A)  £  ♦  i  I 

'  P  r=2 

"  wrS.](q1j^j)Uo-  C  , k )  - 

Hence,  when  we  take  the  Fourier  transform  of  Tld  with  respect  to  w, 
i t  f o 1  lows  that 


3” 1 J  ( k ,  a  ,  , .  .  .  ,  C  n )  :  =  J  nl  dw  T ' J  (  k  ,■  ,w)exp(  -  iw  \)  - 

'  P 

"n  L«V  ''le‘p[ 1  (  kJi  -  J 1 )  -'o  *  '  <  k-  •)  ri  V' J -  i J  (?  O’'  1  ’  r2  ’  ' 


(2.13) 


The  inverse  data  T,d  and  its  analog  TJ 1  are  precisely  what  is  needed 


to  solve  the  inverse  problem  for  q 
,  .  ^  a  i  j  ~  i  i  4] 

a  n  ^  >rr  wii  I  nnr  ^  r  J 


associated  wi th  Eq.  (2.11) 


and  its  analogs  for  ^  ,  see  . 

S,J.  T1^  are  defined  in  terms  ofu^,  respectively .  It  is  poss- 

*-  £  s  l  j 

ible  to  find  a  simple  relationship  between  ^  ,  u  j  which  then  yields 
a  simple  relationship  between  S1^,  T1^.  Actually  if  N  -  3  then  T1^  * 


It  was  shown  in  J  that  the  N-wave  interactions  can  always  be 


reduced  to  three  spatial  dimensions.  It  is  shown  in 
actually  be  reduced  to  two  spatial  dimensions^. 


that  they  can 


RECURSION  OPERATORS  ANO  BI-HAMILTONIAN  FORMULATION  OF  EQUATIONS 
IN  2+1 


.18-20] 


This  is  joint  work  with  P.M.  Santini  Since  it  is  sunmar- 

ized  in  these  proceedings  in  a  separate  contribution,  I  shall  limit 
myself  to  a  few  remarks. 

Since  the  discovery  of  an  exact  approach  to  nonlinear  evolution 
equations  in  1+1  (i.e.,  one  spatial  and  one  temporal  dimension)21^, 
two  interrelated  aspects  have  recieved  much  attention  in  the  litera¬ 


ls  The  development  of  a  method  of  solving  suitable  initial-value 
problems.  For  initial  data  decaying  at  infinity  such  a  method  Is  the 
inverse  scattering  transform  (IST)22^.  This  method  crucially  utilizes 
the  existence  of  an  associated  isospectral  linear  eigenvalue  problem, 
(ii)  The  investigation  of  the  "algebraic"  properties  of  the  given 
equation.  A  fundamental  role  with  respect  to  the  algebraic  properties 
is  played  by  an  integrodifferential  operator,  given  various  names  in 

23 1 

the  literature:  squared  eigenfunction  operator  J,  recursion  opera¬ 
tor24^,  strong  symmetry2^,  hereditary  symmetry2^,  Kahler  operator2*^, 
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regular  operator  J.  This  operator  has  the  following  properties: 

(a)  It  generates  the  associated  hierarchy. 

(b)  It  generates  infinitely  many  symmetries  (in  particular,  if  it  has 
a  certain  property  which  Fuchssteiner  calls  hereditary,  it  gen¬ 
erates  a  set  of  commuting  symmetries) . 

(c)  Its  adjoint  generates  gradients  of  conserved  quantities  (in  part¬ 
icular,  it  generates  a  set  of  involutionary  constants  of  the 


motion) . 

(d)  Its  multiplication  by  one  Hamiltonian  structure  generates  (under 

281 

certain  conditions  J)  a  second  Hamiltonian. 

(e)  The  eigenfunctions  of  its  adjoint  are  quadratic  products  of  the 

eigenfunctions  of  the  associated  isospectral  problem  and  form  a 
29  T 

complete  set  . 

It  should  be  noted  that  given  the  isospectral  eigenvalue  problem, 
there  exists  an  algorithmic  technique  for  obtaining  the  recursion 
operator  (see  for  example^).  This  is,  from  a  unification  point  of 
view,  quite  satisfactory,  since  both  the  method  of  solution  (1ST)  and 
the  algebraic  properties  (recursion  operator)  stem  from  the  same  en¬ 
tity  (isospectral  eigenvalue  problem).  For  the  Korteweg-deVries  (KdV) 

2 

equation  q  =  q  -  6qq  ,  q  =  q(x,t),  the  recursion  operator  0  is  D  - 

u  XXX  X  , 

4q  -  2q  D'1,  where  0*3,  ( D~ Xf ) (x )  */x  f(C)d£.  If  A  is  the  adjoint 
x  x  2  ”  * 

of  ♦,  then  A  satisfies  A'*  -  4Xv  ,  where  g*  solves  <nxx  -  ( q+X ) =  0. 

The  above  two  aspects  have  been  thoroughly  investigated  for  a 
number  of  physically  important  equations  in  1*1.  Each  of  these  equa¬ 
tions  has  physically  significant  two-spatial -dimensional  analogues. 

For  example,  the  KdV  is  generalized  to  the  Kadomtsev-Petviashvi 1 i  (KP) 
equation,  the  modified  KdV  to  the  modified  KP,  the  non-linear  Schrod- 
Inger  to  the  Davey-Stewartson ,  etc.  Furthermore,  these  equations  are 
also  related  to  isospectral  eigenvalue  problems  which  are  appropriate 
generalizations  of  the  corresponding  one-dimensional  ones.  It  is 
therefore  natural  to  investigate  aspects  (i),  (ii)  above  for  two- 
spatial-dimensional  (2*1)  exactly  solvable  equations. 

The  extension  of  the  1ST  to  equations  in  2+1  has  been  recently 
established  in^’6^,  (see  also31^).  However,  the  problem  of  finding 
recursion  operators  in  2+1  has  remained  open;  actually  even  the  exis¬ 
tence  of  such  operators  has  been  doubted  in  the  literature.  In  this 
respect  note: 

1.  The  1ST  of  the  8enjamin-0no  equation  has  all  the  features  of  an 
equation  in  2+1^.  It  is  thus  not  surprising  that  its  recursion  op¬ 
erator  has  not  been  found.  One  of  the  authors  (A.S.F.)  and  8.  Fuch- 
ssteiner,  after  failing  to  find  the  recursion  operator  of  the  Benjamin- 


Ono  equation,  introduced  an  alternative  approach33-^  for  generating 

symmetries.  This  approach  uses  a  certain  function,  called  *;  it  was 

39] 

subsequently  applied  to  a  number  of  equations,  including  the  KP 
However,  for  equations  in  2+1:  (a)  The  relationship  between  t  and  the 
eigenvalue  problem  has  not  been  established,  (b)  There  does  not 
exist  an  algorithmic  way  of  finding  (c)  It  is  not  known  if  -  can 
be  used  to  obtain  the  second  Hamiltonian.  { d )  t  is  not  hereditary. 

2.  The  bi -Hami 1 tonian  nature  of  equations  in  1+1  has  been  emphasized 
as  the  fundamental  property  underlying  i ntegrabi 1 i ty33^ .  However,  the 
bi-Hami 1 tonian  nature  of  all  equations  in  2+1  as  well  as  of  the 
8enjamin-0no  has  remained  open.  The  existence  of  a  recursion  opera¬ 
tor  would  directly  imply  the  second  Hamiltonian,  since  all  these  equa¬ 
tions  have  one  known  Hamiltonian. 

3.  A  number  of  important  results  pertinent  to  the  algebraic  proper¬ 
ties  of  equations  in  2+1  have  obtained  in  the  Soviet  Union36^.  In 
particular  Zakharov  and  Konopel chenko ,  in  a  very  interesting  paper37^, 
claimed  that  recursion  operators  are  purely  one-dimensional  phenomena 

( l . e .  ,  they  do  not  exist  in  more  than  one  dimension).  A  careful  read¬ 
ing  of  their  work  reveals  that  indeed  recursion  operators  of  a  certain 
form  do  not  exist  in  more  than  one  dimension. 

4.  Several  authors  have  noticed  that  mastersymmetries  also  exist 
for  equations  in  1+1;  let  us  call  such  a  mastersymmetry  T.  Actually, 

T  comes  from  a  nongradient  function  and  can  be  used  to  generate  i. 
However,  ;  comes  from  a  gradient  function  and  fails  to  generate  a  re¬ 
cursion  operator. 

The  extension  of  the  inverse  scattering  in  2+1  necessitated  the 
introduction  of  a  new  idea,  the  use  of  a  ( DBAR) .  The  extension  of  the 
theory  of  recursion  operators  and  bi-Hamiltonian  structures  to  equa¬ 
tions  in  2+1  necessitated  the  introduction  of  distributions,  or  more 
precisely  the  introduction  of  integral  representations  of  certain 
differential  operators. 

We  note  that  the  proper  analogue  of  T  is  not  +  but  a  function  de- 
20 1 

noted  in  J  by  T^.  This  function  also  generates  recursion  operators 
in  analogy  with  T. 


4.  THE  AADOMTSEV-PETVIASHVILI  EQUATION  PERIODIC  IN  ONE  SPATIAL  DIM¬ 
ENSION  AND  DECAYING  IN  THE  OTHER 

38l 

This  is  joint  work  with  V.  Papageorgiou  ,  more  details  of  this 
39l 

work  can  be  found  in  .  We  have  considered  four  different  problems: 
API  periodic  in  x,  API  periodic  in  y,  KP 1 1  periodic  in  x,  APII  period¬ 
ic  in  y.  Here  I  will  only  summarize  some  results  for  API  periodic  in 
x  of  period  2 i  . 


4.1.  Analytic  Eigenfunctions. 

The  x-part  of  the  Lax  pair  is  given  by 

iu  ♦  u  +  2i  k  u  =  -u ... 
y  xx  x 

Consider  these  solutions  of  (4.1)  which  also  solve 

r  y 


(4.1) 


t(x,y,k)  =  1+  dn  |  d£  I  +e(m-k,x--;,y-n)u(C  ,n)u'(C,n,k) 

J J.i  mtk *rTL 


-  57  i  dn!  d£  Z  -(m-k,x--;,y--)u(£,n)u"(C,n,k),  (4.2) 

Jy  J-i  mck: :2q 


where 


•/  i-,  ■:  (m.x.y)  =  e 


imx- lm  y 


(4.3) 


These  eigenfunctions  are  9  and  9  with  respect  to  the  complex  k-plane, 
i.e.  they  are  the  boundary  values  of  functions  analytic  in  the  upper 
or  lower  halves  of  the  k-complex  plane. 

A  simple  calculation  shows 


u+(x,y,k)  -  u‘(x,y,k)  -  Z  T(k,m)N(x,y,m,k) ,  (4.4) 

mck+TZ 


where 


T(k,m)  *  57  sgn(k-m)  dn  d£u(  E,  n)u*(C.  n,k)e  ( k-m,  £, n) ,  (4.5) 


X-k,  m-k  are  integer  multiples  of  ;  and  N  solves  the  same  equation  as 
u*(x,y,k)  but  with  the  forcing  replaced  by  8(*-k,x,y). 


4.2.  A  Syranetry  Condition. 

N,  „  are  related  by  (for  simplicity  of  writing  we  suppress  x,y) 
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N(  * ,  k )  =  J  (i  )"(  --k.x.y) 


F(X,m)_  ( m ) 14  (m-k ,  x  ,y )  if  • ->k 


F( > ,m) „"(m)«(m-k ,x ,y )  if  •  ^ k , 


(4.6) 


where  [k+T,*]  =  ;k»T,  k+2x . k+nT-ki,  and 


F(  A,m) 


dn!  dEN(£,n,  k,m)u(E  ,  ") . 
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1  -  i 


(4.7) 


4.3.  Scattering  Equation. 

Using  (4.4)  and  (4.6)  it  follows  that 

u*(k)-u*(k)  *  Z  f(k,m)9(m-k,x,y)u"(m) ;  f  (  k  ,m)=sgn(  k-m)F(  k  ,m) .  (4.8) 
mek+iZ 

Equation  (4.8)  can  be  viewed  as  a  Riemann-Hi 1 bert  (RH)  problem  with  a 
shift  for  u"(k).  The  inverse  problem  consists  of  finding  u*  in  terms 
of  the  scattering  data  f(k,m).  The  time  evolution  of  f(k,m)  is  given 
by 

f (k,m,t)  =  f(k.m,0)e4i(m3'k,)t  (4.9) 

and  f(k,m,0)  can  be  obtained  in  terms  of  initial  data  using  (4.7). 
Remarks 

i)  Localized  solutions,  periodic  in  x  and  decaying  in  y  correspond 

to  homogeneous  solutions  of  (4.2)'.  Such  solutions  have  been 
40 1 

obtained  in  . 

ii)  The  above  formalism  also  follows  by  proper  discretization  of  the 

21 

results  of  J. 

iii)  It  is  interesting  that  a  symmetry  condition  of  the  type  first 

21 

introduced  in  1  is  necessary  not  only  for  KP I  but  also  for  KPII. 

411 

iv)  Some  of  the  above  results  were  first  obtained  by  P.  Caudrey  J 
viewing  KP  as  a  singular  limit  N  *  »  of  a  matrix  N  x  N  one¬ 
dimensional  problem. 

5.  THE  INITIAL  VALUE  PROBLEM  OF  CERTAIN  PAINLEVE  EQUATIONS 

421 

This  is  joint  work  with  u.  Mugan  and  M.J.  Ablowitz  ,  more  de- 

431 

tails  of  this  work  can  be  found  in  J. 


The  mathematical  and  physical  significance  of  the  six  Painleve 
44 1 

transcendents,  PI-PVI  J,  has  been  well  established:  Mathematically, 

these  are  the  only  equations  of  the  form  q  *  F(q  ,q,t),  where  F  is 

rational  in  q^,  algebraic  in  q  and  locally  analytic  in  t  which  have 

the  Painleve  property  (i.e.  their  solutions  are  free  from  movable 

4  5 1 

critical  points).  Physically,  are  closely  related  J  to  physically 
significant  solvable  PDF's  and  have  appeared  in  several  physical 
applications,  see  for  example46"50-^. 

Central  in  the  integrabi 1 ity  of  PDF's  in  1+1  and  2+1  is  their  re¬ 
lation  to  isospectral  eigenvalue  problems.  Similarly,  central  to  the 
integrabil ity  of  the  Painleve  equations  is  their  relation  to  isomono- 
dromic  problems  (see  Sato  et  al51^,  Ueno52^,  Flaschka  and  Newell53-^, 
and  Jimbo  et  a  1 54 ^ ) . 

We  have  systematical ly  considered  the  initial  value  problem  of 

PI  I ,  PIV,  PV.  Equation  PHI  is  contained  in  PV  for  a  special  choice  of 

one  of  the  parameters  of  PV,  equation  PVI  has  been  solved  by  C. 

Cosgrove  and  PI  remains  open.  The  basic  approach  is  that  introduced 
55 1 

in  J  although  we  have  made  certain  simplifications  and  extensions. 
Here  I  briefly  summarize  the  mam  results  using  PIV  as  an  illustrative 

example. 


5.1.  The  Lax  Pair. 


PIV  equation 

d^y  1  ,dy,2  .3  3  ...  2  2  ,  3 

^2  5y  *dT  2  y  4ty  2 ^  u>y  y* 


is  the  compatibility  of  the  following  linear  problems 

t  u  \ 


-(z-9n-9  )  -t 
u  0  00 


f  ( x  ;t) , 


(5.1) 


9 


(5.2a) 


Indeed  Y 


imp] ies 


-4z  ♦  y  ♦  2ty  ♦  4~f 


(5.3a) 


-  2z2 

y 


(Ao  * 


(5.3b) 


£  -  -u(y  *  Zt). 


(5.3c) 


inhere , 


a  =  29 


Eliminating  from  (5.3)  we  obtain  PIV. 

5.2.  Analytic  Eigenfunctions  and  their  Relationship. 

We  recall  that  in  studying  the  initial  value  problem  of  an  equa¬ 
tion  in  1+1  or  in  2+1,  one  uses  the  time-independent  part  of  the  Lax 
pair  to  define  an  inverse  problem,  in  terms  of  certain  scattering  or 
(more  generally)  inverse  data.  Then  one  uses  the  time  dependent  part 
to  find  the  time  evolution  of  these  data.  Similarly,  here  one  uses 
(5.2a)  to  define  an  inverse  probelm  in  terms  of  certain  monodromy  data; 
then  one  uses  (5.2b)  to  fino  the  time  evolution  of  these  data.  To  de¬ 
fine  an  inverse  problem  one  needs  to  consider  the  analyticity  proper¬ 
ties  of  Y ( x ; t )  in  the  whole  complex  x-plane.  Since  Y  satisfies  a 
linear  00E,  its  analyticity  properties  are  completely  determined  from 
the  singular  points  of  (5.2a).  Indeed,  performing  an  analysis  around 
x  =  0,  x=®,  and  introducing  different  solutions  Yj  in  different 

sectors  S.  (so  that  Y.'s  are  normalized  at®)  it  follows  that  (we 
J  J  j 

assume  0  <  9q  <  1 ,  0  <  9^  <  1,  t  ^): 


i)  The  Y^'s,  j  *  1.....5  defined  in  the  sectors  S j ' s ,  where  the 
Sj's  are  given  below  and  each  Sj  contains  the  intial  boundary  line, 
are  related  via 


Yj  +  1  ( * >  ■  Vx)gj*  x  on  CJ  +  1’ 


j  =  1.2.3, 


Y,  ( x )  *  Y4(xe2,7')G,. 


x  on  C, . 


o 


Figure  1. 

6i =  (;  ?)  •  g2 =(s  ?)  ■  g3  -(i  ;)•  G4 = (o  ?)• t5-6 

K  *  G4M»-  =  e'2’"'-). 

1  i )  Yj(x)  Vjx)  *  Y,(x)e0(x)(l/x)°”.  as  !xj-»txinSj,  (5.7 

where  7  (x)  is  nolomorphic  at  x  =  ”, 

2  2 

Q(x)  =  diag[y<-  xt,  -(’j  *  xt)],  =  diagf^,  -ij.  (5.8 

D0 

iii)  Vq(x)  yQ(x)x  ,  as  x  •  0,  x  in  Sj,  (5.9 

where  yg(*)  is  holomorphic  at  x  =  0, 

MQ  =  <jiag(e2  °,  e  °),  y0(xe2’’)  =  yq(x)M0.  (5.10 

iv)  The  connection  matrix  Egis  defined  by 


a0  d0 
y  0  5  0 


yj(x)  =  yq( x ) Eq •  *  in  Sj,  det  Eq=  1. 


From  the  above  it  follows  that  Y  is  a  sectionally  holomorphic 
function.  Its  behavior  at  x  *  0  and  i  =  •  is  determined  from  the  mon 
odromy  matrices  and  M  ;  its  jumps  across  the  Stokes  lines  defined 

in  Figure  1,  are  given  by  the  Stoke's  matrices  Gj . G^.  Hence  its 

entire  behavior  is  determined  from  the  following  data: 

Monodromy  da ta  s  ■  a  ,b,c  ,d ,1q,8q ,  ,  5q  ■ .  (5.12 

5.3.  Properties  of  the  Monodromy  Data. 

It  turns  out  that  all  of  the  monodromy  data  can  be  expressed  in 


iVt’ 


terms  of  two.  In  particular 

V»sEoVEo'  (5'13) 

( l+bc)e^’TTe”  +  [a  ♦  (c>l ) { ab+1 ) ]e_2i7T 6»  =  2cos2tt90.  (5.14) 

Furthermore,  using  (5.2b)  it  follows  that  the  mondromy  data  are  time 
invariant. 


5.4.  The  Inverse  Problem. 

The  inverse  problem  consists  of  reconstructing  Y,  or  more  pre¬ 
cisely  Y  where  Y  *  Ye"^x\  in  terms  of  the  above  monodromy  data. 

The  inverse  problem  can  be  formulated  as  a  Riemann-Hilbert  (RH)  matrix 

problem  along  the  contour  C,  +  C,  +  C,  ♦  C-.  This  RH  problem  is  dis- 
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continuous  both  at  the  origin  and  at  infinity.  Using  the  method  of 
we  separate  this  RH  problem  into  a  sum  of  two  simpler  RH  problems, 
one  defined  on  ♦  C4  and  the  other  on  Cj  ♦  C^.  Furthermore,  it  is 
remarkable  that  the  RH  problem  defined  on  C?  +  can  be  solved  in 
closed  form.  Hence  solving  the  initial  vlaue  problem  of  P1V  is  equi¬ 
valent  to  solving  a  RH  problem  discontinuous  at  x  =  0  and  x  =  «°  and 
defined  on  Cj  ♦  Cj.  This  RH  problem  can  be  mapped  to  a  continuous 
one  using  appropriate  auxiliary  functions.  Hence  its  solution  can  be 
obtained  in  terms  of  a  Fredholm  integral  equation.  Having  obtained  Y 
it  is  straightforward  to  compute  y,  i.e.  the  solution  of  PIV. 


5.5.  Some  Special  Solutions. 

For  certain  choices  of  the  parameters  a,  8,  PIV  admits  one  para¬ 
meter  family  of  solutions  expressible  rationally  in  terms  of  the 
Weber-Hermi te  functions.  Such  solutions  can  also  be  obtained  from  the 
inverse  problem.  Let  9^=  t 9  q  ♦  n  c  J,  then  (5.14),  (5.13)  im¬ 
ply  a  *  c  *  0  and  b  =  -d  respectively .  Let  us  consider  for  concrete¬ 
ness  9  *  -  9  n •  0<9  <  1/2.  Then  one  is  lead  to  the  following  RH 

oo  Q  oo  J 

problem: 

Y +(x)«Y*(x)gy(x) ,  on  Cj+Cj,  9 9 2 1  °n  C3;  g4  on  C 


(5.15) 


r 


*mmm$m**i*  m  mi*  imp««u  m  w  m  tmmmmimmm**#* 


This  RH  problem  can  be  solved  in  closed  form: 


•  -0  e 

®  CO 

x  x  F(x 


d  !  e^(x)  „  2 

;  F(x)  *  jr-j-  j  x  2  *  — - dx,  q(  x )  =  -y  *  xt, 

'  L  x-x  c 


(5.16) 


where  the  contour  L  is  given 


Note  that 

F(x)  *  -  ^f(t)  +  0(y),  f(t)  =  He^-  _  ^  ( i  t ) ,  He  *  Weber-Hermi te 
’’  x  i" 

function. 

5.6.  Schlesinger  Transformations . 

These  transformations  change  9Q,  9^  to  9q,9J,  where 


a  1  =  p  ♦!!!  a'rG 

-  ?*  '3 


i,  m,  n  c  Z. 


(5.17) 


To  obtain  these  transformations  let  Y  *  R(x,t)Y  correspond  to  9q,  9^ 
but  to  the  same  monodromy  data.  Then  it  can  be  shown  that  R  satisfies 
*  very  simple  RH  problem,  which  can  actually  be  sovled  in  closed  form: 
for  example 


i 


90  ’  90  *  2  ' 


9‘  *  9  ♦  i  , 
»  2 


y(z*V- 


(5.18) 


6.  SOLUTIONS  TO  A  CLASS  OF  NONLINEAR  SINGULAR  INTEGRO-OIFFERENTIAL 
EQUATIONS 

This  is  joint  work  with  M.J.  Ablowit*  and  M.D.  Kruskal^  and  was 


motivated  by  some  recent  results  of  Constantin,  Lax,  and  Majda  , 


these  authors  used  the  transformation  w 


u  ♦  iHu,  where  H  denotes  the 

lj..  ,  ~  r\r\r  _  1  2 


Hilbert  transform,  to  map  the  equation  ut  =  uHu  to  the  ODE  w  =  -  jw . 

Using  the  fact  that  w  =  u  *  iHu  is  the  boundary  value  of  a  func¬ 
tion  analytic  in  the  lower  half  complex  x-plane  it  follows  that 
Hw  -  -iw,  and  more  generally  Hwn  =  -iw11,  Hew  =  -iew  etc.  This  enables 
us  to  map  a  large  class  of  nonlinear  singular  i ntegro-di f ferent ia 1 
equations,  via  explicit  transformations ,  to  either  ordinary  different¬ 
ial  equations  or  to  linearizable  partial  differential  equations.  Con¬ 
versely,  given  a  linearizable  PQE  there  is  an  algorithmic  way  of  find¬ 
ing  its  singular  integro-di fferential  analogue.  Examples  include: 


(a)  ODE’S 


(b)  Singular  Inteqro-Di fferential  Equations 


u^  =  uHu 


ut  -  u 


3u(Hu)‘ 


Hu 

u ^  '  e  sm  u. 


(a)  p0E's  (b)  Singular  INtegro-Di fferential  Equations 

wt  =  wxx‘i(w2)x  ut  °  uxx  +  ? ( uHu ) x  (6.4) 

VWxxx-ia(w2)x*g(w3)x=°  ut+uxxx*2a(uHu)x  +  3(u3-3u(Hu)2)x  =  0  (6.5) 

«xt  =  ie"1w  uxt  =  eHusin  u  (6.6) 

wt  *  '("xx+("2)x)=0  ut  '  ^Hu^xx  +  2(uHu)x-  (6.7) 

t^-(w +w  -i(w^)  )=-3c^w  r“ ( u ,  +u  +(2uHu)  )  =  -3o^u  .  (6.8) 

ax'  t  xxx  '  yy  9x'  t  xxx  '  %  yy  '  ' 

Equation  (6.4a)  is  essentially  8urgers  equation  and  can  be  linearized 
via  the  Cole-Hopf  transformation  w  =  - i ( inf ) x .  Equation  (6.4b) 
arises  in  various  population  ecological  models  and  to  our  knowledge, 
was  first  considered  and  solved  via  a  dependent  variable  transforma- 

CO  "I 

tion  and  splitting  into  upper  and  lower  functions  by  J.  Satsuma  . 

In  equations  (6.5)  a,  6  are  real  constants,  and  (6.5b)  is  an  analog 
of  the  Gardner  equation  (a  combination  of  KdV  and  modified  KdV). 
Equation  (6.6b)  is  related  to  the  Liouville  equation  (6.6a)  and  is 
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known  to  be  1 i nearizable . 

A  (3+l}-dimensional  equation  can  also  be  linearized  via  (6.8a). 
Namely  let  H^u  denote  the  Hilbert  transform  of  u(x,y,z,t)  with  respect 
to  the  variable  z,  i .e. , 

h  u .  i  r  dc. 

2  ’T  !  -2 

'  —  3C 

Then  instead  of  K-P:  we  may  consider  a  mul ti -dimens lona 1  analog  of 
(6.8b) 

~-(u„  ♦  u  +  2 ( uH  u)  )  =  -3c^u  ,  (6.9) 

3x'  t  xxx  '  z  'x'  yy 

and  it  is  also  mapped  to  the  KP  equation  (6.8a),  via  w  =  u-mH^u. 

7.  HODOGRAPH  TRANSFORMATIONS  OF  LINEARIZABLE  POE'S 

59  T 

This  is  joint  work  with  P.  Clarkson  and  M.J.  Ablowitz  .  Since 
P.  Clarkson  will  present  these  results  in  a  separate  contribution  of 
these  proceedings,  I  will  only  make  a  few  comments. 

We  call  two  PDF's  equivalent  if  one  can  be  obtained  from  the 
other  by  a  transformation  involving  the  dependent  variables  u  *  :  (v) 
and/or  the  introduction  of  a  potential  variable  (u  =  vx  or  ux  =  v). 

It  is  well  known60 ^  that  the  most  general  semi  1  inear  equation  of 
the  form 

ut  ‘  uxx  *  f(u’ux)  {7A) 

which  is  1 inearizable,  is  either  linear  or  the  Burger's  equation 
(which  with  the  above  definition,  is  equivalent  to  a  linear  equation). 
Fokas  and  Yortsos61^  have  shown  that  the  most  general  quasi  1  inear  equa¬ 
tion  of  the  form 

ut  -  g(u)uxx  ♦  f(u,ux)  (7.2) 

which  is  1 inearizable,  is  equivalent  to  the  equation 

u{  *  (u  ux)x  +  au  ux,  a  is  an  arbitrary  constant.  (7.3) 

The  above  equation  can  be  mapped  to  the  Burger’s  equation  via  an 
extended  Hodograph  transformation,  i.e.  a  transformation  of  the  form 

t  =  t,  C  *  j  u ( x * ,t)dx* .  (7.4) 


It  Is  also  well  known  that  the  Harry-Dym01- J  equation,  can  be  mapped 
to  a  modified  Korteweg-deVries  (MKdV)  equation  via  an  extended  Hogo- 
graph*^.  In  that  sense,  the  Harry-Dym  equation  is  a  quasilinear 
analogue  of  the  MKdV  equation. 

One  is  naturally  motivated  to  ask  the  following  questions: 
i)  Is  there  an  algorithmic  method  of  finding  a  quasilinear  analogue 
of  any  semil inear  equation? 
i  i )  Is  the  associated  quasilinear  equation  unique? 

i l i )  Conversely,  given  a  quasilinear  equation,  is  there  an  algorithmic 
method  of  finding  whether  it  can  be  mapped  to  a  semilinear  equa¬ 
tion  as  well  as  finding  this  semilinear  equation? 

In  the  above  paper  we  consider  the  above  questions  for  semilinear 
and  quasilinear  equations  of  the  type 

ut  =  Unx  *  f(u’ux . u(n-l)x}>  n  -  2’  Unx  t  ~n  '  (7-5) 

and 

ut  =  9(«)unx*f(u,ux . Vl)*)’  n  -  2'  S  ’  °  (7‘6) 

respectively.  The  answer  to  question  i)  above  is  affirmative.  Also, 
the  associated  quasilinear  equation  is  unique,  in  the  sense  that  ex¬ 
tended  and  pure  hodograph  transformations  yield  equivalent  quasilinear 
equations ■  (By  pure  hodograph  we  mean  transformations  of  the  form 
-  -  t,  £  -  u(x,t)).  Furthermore,  we  find  the  most  general  equation 
of  the  form  (7.6)  which  can  be  mapped  via  an  extended  hodograph  trans¬ 
formation  to  a  semilinear  form. 

The  above  results  are  of  some  interest  in  establishing  whether  an 
equation  is  a  candidate  for  linearization.  Suppose  that  one  is  inter¬ 
ested  in  investigating  whether  a  given  quasilinear  equation  is  linear- 
izable.  We  propose  the  following  algorithmic  procedure: 

1.  Put  the  equation  into  its  potential  canonical  form 


V^V  +  H(v ,,V  .. ■ • .V 

x  nx  x  xx 


by  using  the  transformation  v^  =  g 


(n-l)x 

-1/n 


), 


(7.7) 


(u). 


Apply  a  pure  hodograph  transformation  to  equation  (7.7).  If 
equation  (7.7)  is  transformable  to  a  semilinear  equation,  it  will 
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become 

ut  =  unx  *  F(W . u(n-l)x) '  (7'8> 

3.  Investigate  whether  equation(7.8)  is  1 i neari zabl e .  This  is  eas¬ 
ier  than  investigating  whether  (7.6)  is  linearizable  directly. 

The  reason  for  this  is  twofold.  First,  for  at  least  third  order 
equations  there  is  a  complete  classification  of  all  linearizable 
equations.  Within  equivalence,  there  exist  only  six  such  equa- 
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tions  .  Hence  one  simply  needs  to  study  if  there  exists  an 
equivalence  transformation  to  map  equation  (7.8)  with  n  =  3,  to 
one  of  these  six  canonical  equations.  Second,  for  equations  with 
n  >  4  one  may  investigate  the  question  of  linearization  via  the 
Painlev6  test.  We  point  out  that  quasilinear  partial  different¬ 
ial  equations  do  not  appear  suitable  for  applying  the  Painlevg 
test.  Ramani,  Dorizzi  and  Grammaticos65^  (see  also67^  and  the 
references  therein),  introduced  the  notion  of  "weak-Pa i nlevg"  in 
order  to  deal  with  equations  such  as  the  Harry-Dym  equation  which 
are  linearizable  after  a  change  of  variables.  However,  the 
higher  KdV  equation 


u»  =  u 
t  xxx 


3 

u  u 


although  it  is  thought  to  be  nonl inearizable  (it  has  only  three 
independent  polynomial  conservation  laws  of  a  certain  type66^) 
and  is  also  "weak-Painlev£"®7^ .  Therefore  the  "weak-Painleve" 
concept  does  not  seem  to  distinguish  between  a  linearizable  and 
a  nonl inearizable  equation. 

We  point  out  that  one  often  finds  in  the  literature  claims  of 
"new"  third  order  linearizable  equations.  These  equations,  using  the 
notion  of  equivalence  can  be  mapped  via  a  pure  hodograph  transforma¬ 
tion  to  one  of  the  six  canonical  equations  mentioned  above. 

The  above  algorithmic  approach  is  useful  provided  that  a  given 
linearizable  quasilinear  equation  can  be  mapped  to  a  semilinear  form. 
The  above  approach  will  fail  if  there  exist  linearizable  quasilinear 
equations  which  can  not  be  mapped  to  a  semilinear  form.  It  is  shown 
in^  that  such  equations  do  not  exist  for  at  least  n  -  2.  The  ques- 
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tion  of  whether  such  equations  exist  for  n  _>  3  remains  open. 

8.  THE  SCALING  REDUCTION  OF  THE  THREE-WAVE  RESONANT  SYSTEM  AND  THE 
PAINLEVE  VI  EQUATION 

This  is  joint  work  with  R.A.  Leo,  L.  Martina,  and  G.  Soliani8®-^. 
It  is  well  known  that  for  a  large  class  of  equations,  the  large 
time  asymptotic  limit  is  governed  by  certain  similarity  solutions  of 
the  underlying  POE.  If  this  PDE  is  an  exactly  solvable  equation  in 
1  +  1  (i.e.  in  one  spatial  and  in  one  temporal  dimension)  one  expects^ 
that  the  similarity  solutions  satisfy  an  ODE  of  the  Painlevfc  type. 

In  the  above  paper  we  considered  the  three-wave  resonant  inter¬ 
actions  in  the  case  of  explosive  instabi  1  i ty6^*2^ . 

ujt  +  cjujx  ’  lu*uk  =  °’  1’2'3*  >>  k  *1'  (8'1) 

where  u^-fx.t)  are  the  complex  amplitudes  of  the  wave  parameters,  c^ 
are  their  group  relations  and  *  denotes  complex  conjugate.  Assuming 
C1  <  c2  <  c3  and  us’n9  invariance  under  x-translation,  t-translation 
and  appropriate  scaling  we  are  lead  to  consider  a  system  of  three 
first  order  nonlinear  ODE'S. 

This  system  via  a  series  of  transformations,  can  be  mapped  to  a 

single  second  order  ODE,  which  is  quadratic  in  the  second  derivative. 

Such  equations  are  outside  the  class  investigated  by  Painlevg  and  his 
44 1 

school  J,  however  the  equation  obtained  above  is  a  particular  case  of 
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an  equation  recently  studied  by  Fokas  and  Yortsos'  ,  in  their  invest¬ 
igation  of  exact  trasnformations  of  Painleve  VI  equation: 

A  fundamental  role  in  the  exact  treatment  of  the  Painlev^  equa¬ 
tions  is  played  by  certain  transformations  which  map  solutions  of  a 
given  Painlev£  to  solutions  of  the  same  Painlev£  but  with  different 
choice  of  the  parameters.  Such  transformations  for  PII-PVI  were  given 
in2*  respectively.  Finding  such  a  transformation  for  PVI  necessi¬ 
tated  the  introduction  of  an  auxiliary  equation  which  is  quadratic  in 
the  second  derivative.  However,  it  was  shown  in7^  that  this  equation 
can  be  mapped  to  PVI . 

The  second  order  ODE  obtained  from  the  similarity  reduction  of 
the  three-wave  resonant  interactions  is  a  special  case  of  the  above 


auxiliary  equation  (see  equation  (2)  of 
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